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Abstract . We present the basic princi- 
ples of very- long-baseline interferometry 
as related to its use in the determina- 
tion of vector baselines, polar motion, 
and earth rotation. 

Introduction 

Ten years ago, almost to the day, the 
first successful bandwidth-synthesis VLBI 
measurements were made. It is thus appro- 
priate to now review the principles under- 
lying the technique. The review will be 
restricted to aspects relevant to geodet- 
ic applications that involve observations 
of extragalactic radio sources. In such 
applications, arrays of two or more radio 
telescopes observe any given source simul- 
taneously. From sets of observations of 
a suite of such sources, one can obtain 
the desired geodetic information: base- 
line-vector, polar-motion, and earth- 
rotation parameters.* 

We shall first describe briefly the in- 
strumentation used in these observations 
and then discuss the basic observables 
and their simplest interpretation. Fi- 
nally, we consider some complications 
of the interpretation due to the various 
geophysical "signals" and non-geophysical 
"noise" that affect the observables. 

Instrumentation 

A VLBI system consists of an array of 
at least two antennas that observe the 
same radio source simultaneously. A di- 
rect electrical connection is not main- 
tained between the antennas, thus allow- 
ing them to be separated by thousands of 
kilometers. The local-oscillator signals, 
used at each antenna to convert the radio- 
frequency signals from the source to the 
video (low-frequency) band, are derived 
from a frequency standard at the site. 
These standards are sufficiently stable 
that the relative phases of the signals 
from the source received at the two an- 
tennas are preserved. 

*Several groups are currently engaged in 
such geodetic applications of VLBI: (1) 

A group from the Goddard Space Flight Cen- 
ter, the Haystack Observatory, the Massa- 
chusetts Institute of Technology, and the 
National Geodetic Survey; (2) a Jet Pro- 
pulsion Laboratory - National Geodetic 
Survey group; and (3) a Canadian-British 
collaboration. A European consortium is 
also being organized for similar purposes. 

Proc. of the 9th GEOP Conference. An huermuional Sympimum on the AppUcmUms of 
Geodesy loGeodyiwmics. October 2-5. 1978. Dept, of Geodetic Science Rept. No. 280. The 
Ohio State Univ.. Columbus, Ohio 43210. 


The video signals are recorded on mag- 
netic tape at each site, with the refer- 
ence time for the recordings being de- 
rived from the same standard as is used 
to govern the local-oscillator signals. 

The tape recordings are then transported 
to a common center where those recorded 
simultaneously are cross-correlated to ob- 
tain the basic VLBI observables. 

Observables 

The basic observables in geodetic VLBI 
experiments are (i) the difference in the 
times of arrival at two antennas of a sig- 
nal from a source; and (ii) the rate of 
change of this time difference. The 
measurement of the time-of-arrival dif- 
ference can be of two types : the phase- 
delay difference or the group-delay dif- 
ference. The phase-delay-difference ob- 
servable (hereinafter "phase delay"), 
being based on measurements of phase, can 
be obtained very precisely, but usually 
ambiguously, due to the inability to re- 
solve the "2TTn" problem. The group-delay- 
difference observable (hereinafter "group 
delay") is usually determined with less 
precision than is the phase delay. But 
the group delay, being determined by the 
rate of change of phase delay with fre- 
quency, is usually unambiguous. To ob- 
tain reasonable accuracy in the measure- 
ment of group delay it is necessary to 
make phase-delay measurements over a wide 
band of frequencies simultaneously, or 
nearly simultaneously, since the uncer- 
tainty in the group-delay measurement is 
inversely proportional to this bandwidth. 
In practice, only a relatively narrow 
band of frequencies can be recorded. But 
this band can be split up into narrower 
bands which are spread over a very wide 
band. This technique is called bandwidth 
synthesis. In the newest, Mark III, sys- 
tem (see below), 28 narrow bands, each 2 
MHz wide, are being distributed over a to- 
tal band of up to 400 MHz. If the error 
in the measurement of phase for any one 
band is a(ip), then the error a(Tg) in the 
measurement of the group delay will be 
given by o(Tg) ~ a(ip)/Af, where Af is the 
rms spread of- the center frequencies of 
the individual bands about their mean. 
Because these individual bands do not cov- 
er the entire spanned band, the estimate 
of the group delay, too, could be ambigu- 
ous. However, a proper choice of the 
spacing of the individual bands, as ex- 
plained below, can insure that any inher- 
ent ambiguity can be eliminated reliably. 

The difference of the phases of the 
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signals that would be received at two an- 
tennas as a function of frequency exhi- 
bits curvature due primarily to the ef- 
fects of the earth's ionosphere: the low- 
er the frequency, the shorter the phase 
delay. The actual shape of the phase vs. 
frequency curve for VLSI observations 
will depend on the relative amounts of 
plasma between the source and the two 
sites. But with appropriate spacings of 
the narrow bands, one can "connect" un- 
ambiguously the phase at one band with 
those at all the other bands. For exam- 
ple, bands spaced in accord with a geo- 
metric series allows use of a bootstrap 
technique to connect phases first between 
the closest bands and then between the 
more widely separated bands through use 
of the characteristics of the curve es- 
tablished by the connection between the 
closest bands. Of course, the "absolute" 
phase would still be uncertain by mul- 
tiples of 2 tt, but the relative phases be- 
tween bands would be freed from any such 
ambiguity. Only the relative phases af- 
fect the group delay which is equal to 
the slope of the curve of phase delay vs. 
(angular) frequency. 

Information Content 

We now consider the information con- 
tent of the observables under simplified 
assumptions. In particular, let us ig- 
nore the propagation medium and assume 
that the earth is rigid and rotates with 
a constant, known, angular velocity. We 
may then write the expression for the de- 
lay observable as a function of time as: 

x(t) = il(t)«i + [t-to] , (1) 

where B is the baseline vector connecting 
a pair of antennas, s is a unit vector in 
^he direction of the source, and and 
x^^ are, respectively, the offsets in 
epoch and rate of the clock at one site 
with respect to those at the other site. 
This equation represents a diurnal sinu- 
soid added to a straight line. The first 
term in the equation contributes the di- 
urnal sinusoid due to the rotation of the 
baseline vector in inertial space. The 
slope of the straight line is due to the 
clock-rate offset and the intercept is 
due to a combination of the clock-epoch 
offset and j^he product of the polar com- 
ponents of B and s. Clearly this curve 
can be specified by four parameters: 
the intercept and slope of the straight 
line, and the amplitude and phase of the 
sinusoid. (We do' not consider the period 
of the sinusoid, since that is given by 
assumption. ) Thus four measurements of 
the delay suffice, in principle, to de- 
termine x(t); any additional measurements 
will be redundant. But how many unknown 
parameters are there for this situation? 


Naively, one would conclude that the base- 
line vector contributes three, the source 
two, and the clocks two, for a total of 
seven. However, the origin of right as- 
cension of our system is arbitrary; only 
the origin of declination is fixed by the 
assumption of a known angular velocity 
for the earth. Since the right ascension 
of the source can be used to define this 
arbitrary origin, the number of unknown 
parameters is only six. Nonetheless, a 
unique solution cannot be obtained for 
these six parameters from observations of 
a single source. Observations of each ad- 
ditional source adds two unknowns; the 
coordinates of the source on the plane 
of the sky. But such observations 
also can be used to determine three ad- 
ditional parameters: the intercept of the 
straight line and the amplitude and phase 
of the sinusoid appropriate for the ad- 
ditional source. The slope of the 
straight line provides no new information 
since it is determined solely by the 
clock-rate offset. It is clear that with 
four observations of one source and three 
each of two more sources, a useful solu- 
tion for all of the parameters of this 
simple model can in general be determined. 
(The measurements of delay rates simul- 
taneously do not reduce the requirement 
for observations of three sources.) The' 
accuracy of the determination of these 
ten parameters will depend, of course, 
not only on the accuracy of the measure- 
ments of delay (and delay rate) , but also 
on the baseline, the distribution of the 
sources in the sky, and the distribution 
of the observations in time. 

Complications 

The situation actually encountered 
with VLBI is, of course, far more compli- 
cated than outlined in the previous sec- 
tion. We can conveniently divide these 
complications into two categories: sig- 
nals and noise. Here signals refer to 
those effects on the observables which 
are of geophysical interest, and noise 
refers to those of no intrinsic interest. 
(This point of view, needless to say, is 
a rather parochial one since one person's 
noise is often another's livelihood.) 

We shall consider precession, nutation, 
solid-earth tides, crustal motions, varia- 
tions in UTl, and polar motion to be sig- 
nals. On the other hand, clock instabil- 
ities and uncertainties in our knowledge 
of source characteristics and of the pro- 
pagation medium shall be considered as 
noise. We discuss each set in turn. 

Signals 

Precession and Nutation . Changes in 
the direction in space of the spin axis 
of the earth with periods long compared 
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to a day are sensed with VLBI through the 
corresponding changes in the coordinates 
of the radio sources. These changes will, 
however, preserve the arclengths between 
sources. At present, estimates of the 
precession constant made from VLBI meas- 
urements have an uncertainty of a few 
tenths of an arcsecond per century , several- 
fold larger than the uncertainty associ- 
ated with the presently accepted value 
based on optical observations. The VLBI 
estimate is consistent with the optical 
one to within twice the formal standard 
error of the former. No estimates of any 
of the nutation terms have yet been made. 

Solid-Earth Tides . The semi-diurnal 
solid-earth tide imparts a distinctive 
signature to the VLBI observable, since 
almost all other effects introduce a di- 
urnal signature. The maximum magnitude 
of this effect on the observable has been 
slightly greater than one nanosecond. 

Thus estimates of the local values of the 
vertical and horizontal Love numbers, 
and h, can be obtained from VLBI data. 
Current estimates agree, to within their 
uncertainty of cibout 0.05, with the "ex- 
pected" values. 

Crustal Motion . Changes in crustal 
configuration can be sensed by long-term 
changes in baseline lengths; in addition, 
significant changes in the corresponding 
baseline directions for an array of an- 
tennas would signify crustal motions pro- 
vided that these changes were incompati- 
ble with a rigid rotation of the array. 

No measurements of crustal motions have 
yet been obtained from VLBI data, but 
there is every reason to believe that 
such motions will be detected within a 
few years. 

UTl and Polar Motion . Variations in 
the rate of rotation of the earth and in 
the position of the axis of figure with 
respect to the axis of rotation affect 
the directions of baseline vectors. As 
with the arclengths between sources with 
respect to precession and nutation, the 
lengths of baselines are unaffected by 
such variations in the rate of rotation of 
the earth and in the position of the pole. 

We should stress that the VLBI observ- 
ables, for an arbitrary baseline, have no 
sensitivity to the "initial" orientation 
of the earth and direction (in space) of 
its axis of figure. Only changes in these 
quantities can be detected. In addition, 
any "common-mode" errors in the epoch set- 
tings of the clocks at the antenna sites 
will be indistinguishable from correspon- 
ding changes in the orientation of the 
earth about its spin axis (UTl) . Finally, 
note that VLBI data obtained for one base- 
line are sensitive to only two indepen- 
dent combinations of the three parameters 
needed to specify changes in the position 
(in space) of the axis of figure of the 
earth and in the orientation of the earth 


about this axis: Such changes affect 

only the direction of the baseline which 
is described by only two independent 
parameters. Consider, as an example, a 
wholly north-south baseline. With such 
a configuration, the VLBI data would have 
no sensitivity to a rotation of the earth 
about the pole. For an east-west base- 
line, on the other hand, polar motion in 
a direction along the meridian of the mid- 
point of the baseline would not affect 
the VLBI observables. In both these 
cases, the baselines would undergo paral- 
lel displacements which cannot be detect- 
ed from observations of sources "at in- 
finity." Two baselines, or at least three 
antennas, are needed ,in a VLBI array to 
detect all three components of the chan- 
ges in UTl and pole position. These base- 
lines must, of course, not be parallel. 

Estimates of UTl and polar motion from 
VLBI data now have accuracies comparable 
to those of other techniques, such as the 
classical optical methods, the Doppler 
tracking of satellites (for polar motion) , 
and the laser ranging to retrof lectors on 
the moon (for UTl) . It is expected that 
the accuracy of the VLBI estimates will 
improve nearly tenfold within the next 
five years due primarily to the intro- 
duction and use of the new, Mark III, 

VLBI system. A prototype of this system 
has already been tested successfully; 
five copies of the complete system are 
currently under construction for place- 
ment at suitably distributed antennas. 

It will be important to check the im- 
proved determinations of polar motion and 
UTl through redundancy and through com- 
parison with the results from other im- 
proved techniques in order to assess the 
accuracy of these determinations. 

Noise 

Clock Instabilities . The two parame- 
ters, for clock epoch and rate offsets, 
do not provide an adequate representation 
of the relative behavior of the clocks at 
any two antennas of an interferometric 
array over the period of many hours need- 
ed to determine the baseline vector, 
source positions, etc. This statement 
applies to the current field units of all 
atomic clocks, including the hydrogen- 
maser frequency standards. 

A number of possibilities exists to 
minimize the impact of these clock insta- 
bilities on the accuracy with which geo- 
physical information can be extracted 
from VLBI data. First, one can use high- 
er-order polynomials to represent the 
relative clock behavior; here the point 
of "diminishing returns" sets in at about 
the sixth order. Second, the clock per- 
formance, especially of hydrogen-maser 
standards , can be improved to match that 
achieved in the laboratory. The sensiti- 
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vity of the maser standards to environ- 
mental effects can also be reduced to 
minimize the introduction of diurnal sig- 
natures into the VLBI data. Third, one 
can reduce the effects of long-term 
drifts in the relative clock behavior by 
using clock stars. Thus, one can make 
observations repeatedly, say every hour, 
of some suitable source and use these ob- 
servations to correct for the relative 
clock drifts. To be suitable, this 
source should be visible from both sites 
for a large fraction of the diurnal cy- 
cle and should yield a reasonably large 
correlated flux density so that accurate 
delay observations are possible to make. 

Source Characteristics . The radio 
sources affect the determination of geo- 
physically interesting quantities through 
the strength of their radio emissions, 
their distribution on the sky, and the 
accuracy with which we can determine 
their positions. These positions, in 
turn, depend on the structure and inter- 
nal kinematics of the regions of radio 
emission in each object. 

At present, the entries in the catalog 
of known, and potentially-usable , extraga- 
lactic radio sources number in the hun- 
dreds. Positions of a few dozen of those 
with the strongest emissions are now be- 
ing determined routinely with an estima- 
ted accuracy of about 0V02, except for 
the declinations of sources that lie near 
the equatorial plane. The accurate de- 
termination of the declination of those 
sources requires the use of interferome- 
ters with baselines that possess large 
components in the north-south direction. 
Few such baselines have so* far been avail- 
able for extensive sets of measurements. 

Aside from the examination of the char- 
acteristics of the postfit residuals, the 
main method for assessment of the accura- 
cy of source-position determinations is 
the comparison of results obtained with 
different equipment. Such comparisons, 
made several years ago and based on data 
obtained with somewhat less advanced VLBI 
systems, showed agreement to within cibout 
O’.'OS rms. (Note that a O'.'OOl error in 
source position corresponds approximately 
to a two-centimeter error in length for a 
4,000-km baseline.) 

Most extragalactic radio sources are 
not "points" when viewed on the scale of 
milliarcseconds . Rather, they exhibit 
complicated structure. This structure 
in their brightness can be mapped and a 
suitable feature in the map, or the over- 
all ■ center. of brightness , can be used as 
a reference point. There-are, however, 
technical- difficulties: in: the determi- 
nation of: unambiguous brightness maps. 
These dif f icul-ties -are being overcome 
and reliable maps on the scale of tenths 
of a milliarcsecond are now being ob- 
tained for some of -the radio sources. 


There is yet a further difficulty in 
the use of extragalactic radio sources: 
most are not static. Dramatic changes 
have been observed in the brightness 
structure of some of these sources at 
the level of tenths of a milliarcsecond 
in angular resolution on a time scale of 
a few months. Thus, to enable positions 
of extragalactic radio sources to be used 
effectively as a reference system at the 
level of milliarcsecond accuracy for geo- 
physical applications of VLBI, one must 
monitor the brightness distributions of 
these sources as a function of time and, 
perhaps, as a function of radio frequency 
as well. 

Propagation Medium . In regard to any 
substantial effects on VLBI data, the pro- 
pagation medium can be considered to be 
composed of two components : the iono- 
sphere and the troposphere. The effects 
of the ionosphere can, and will, be virtu- 
ally eliminated by observing simultaneous- 
ly in two widely separated radio frequen- 
cy bands (r2GHz and rSGHz) , The Mark III 
VLBI system is equipped for such dual- 
band observations. Moreover, enough suit- 
able sources exist to allow effective use 
of the dual-band technique. 

The troposphere is in effect non-dis- 
persive at radio frequencies and is there- 
fore a more troublesome contributor of 
noise. The troposphere can also be de- 
composed into two components: wet and dry. 
For the latter, the assumption of hydro- 
static equilibrium is a very good one; 
measurements at each site of surface pres- 
sure combined with a good model of the at- 
mosphere, then allows a good estimate to 
be made of the phase delay added in the 
zenith direction by the dry component: 
about 7.5 nanoseconds (equivalent to an in- 
crease in path length of about 2.3 m). It 
is widely thought that the error in this 
estimate can be kept at the 0.1% level or 
perhaps below. Mapping to other zenith 
angles, however, will increase the error 
somewhat since the "slant" . atmospheric 
path length cannot be determined so ac- 
curately from the pressure measurement 
at -the antenna site. The situation with 
the wet component is more difficult. The 
water-vapor in -the atmosphere is not in 
hydrostatic equilibrium and is quite 
variable in amount. Although the total 
effect on -the path, length of radio waves 
is, on average, only about 7% of that of 
the dry component, the wet component can- 
not be modeled accurately. Various sim- 
ple techniques have been used to try to 
ameliorate this problem. Such techniques 
involve various combinations of model at- 
mospheres and mapping functions with or 
withou-t dependence on surface measure- 
ments of temperature, pressure, and dew 
point, and with or without parameters 
that can be estimated for each site from 
short, ~8 hr, spans of data. Unfortunate- 
ly, these techniques may well be defi- 
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cient, especially for long baselines, in 
removing the effects of the atmosphere 
on the estimates of the "vertical" com- 
ponent of the baseline with which they 
are highly correlated. 

The technique which has elicited the 
greatest expectations for providing the 
solution to the wet-component problem is 
based on the use of radiometer measure- 
ments at each site of the brightness tem- 
perature of the atmosphere at and near 
the =^23GHz line in the spectrum of water- 
vapor emission. Studies indicate that 
this brightness temperature can be rela- 
ted with reasonably high accuracy to the 
excess path length attributable to the 
water-vapor content of the atmosphere. 
However, to date, almost all VLBI results 
have been obtained without the benefit 
of water-vapor radiometer measurements. 

Atmospheric effects thus loom as the 
limiting factor in the accuracy achiev- 
able with VLBI in the determination of 
geophysical quantities. What will that 
limit be? An assessment based on theory 
alone is unlikely to be accurate. Meas- 
urements are clearly called for. Series 
of VLBI experiments should be made with 
supplementary water-vapor radiometer 
measurements, under a variety of local 
weather conditions, and for various base- 
line lengths. For short baselines, up 
to several kilometers in length, inde- 
pendent determination of the baseline 
vector can usually be made, with some 
effort, at the millimeter level of ac- 
curacy by means of conventional survey 
techniques. For long baselines, up to 
several thousand kilometers in length, 
independent means of verification at the 
relevant level of accuracy seem to be 
limited to laser ranging to artificial 
satellites or to the moon; such verifi- 
cation, however, will not be easy nor 
inexpensive for a number of practical 
reasons. The repeatability and consis- 
tency of VLBI results themselves may 
well have to provide the main standards. 
Since suitably accurate, and independent, 
estimates of UTl and polar motion may 
not be available, repeated checks on the 
individual components of the vector base- 
line will likely require multi-site ex- 
periments, say with four or more separa- 
ted antennas, to reduce the confusion be- 
tween UTl and pole position changes on 
the one hand and changes in baseline di- 
rection on the other. The many antenna 
sites serve to over ^determine UTl and 
polar motion, with the redundancy pro- 
viding the meaningful check on the con- 
sistency of the estimates of some of the 
baseline components. For some combina- 
tion of the precision and the time 
spanned by the sets of measurements, one 
must be concerned also about the genuine 
changes in baselines expected from plate 
tectonics; of course, detection of such 
changes are a major purpose of the meas- 


urements. 

At present, checks on the repeatabili- 
ty of baseline determinations have in- 
volved primarily two-element interfero- 
meters. For short baselines, of the or- 
der of one kilometer in length, repeata- 
bility has been obtained at the five mil- 
limeter level in all components, and 
verified later by the results of a con- 
ventional survey. For long baselines, 
of the order of several thousand kilo- 
meters in length, only repeatability in 
baseline length has been meaningful; here 
the spread about the mean in a recent 
series of a dozen sets of measurements 
was under five centimeters. The source 
positions used in the analysis of each 
of these experiments were fixed in ac- 
cord with the results from the ensemble 
of experiments; errors in these positions 
tended therefore not to have a serious 
effect on the repeatability of the de- 
terminations of baseline length. 

Conclusion 

The future of VLBI as applied to geo- 
detic and geophysical problems, especial- 
ly to the determination of UTl and polar 
motion, looks quite bright. Although the 
last ten years have been devoted almost 
exclusively to the development of VLBI, 
the next ten should yield significant 
results. 
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